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Abstract 

Generalized VT symmetry provides crucial insight into the sign problem for two classes of models. 
In the case of quantum statistical models at non-zero chemical potential, the free energy density 
is directly related to the ground state energy of a non-Hermitian, but generalized PT-symmetric 
Hamiltonian. There is a corresponding class of "PT-symmetric classical statistical mechanics models 
with non-Hermitian transfer matrices. For both quantum and classical models, the class of models 
with generalized VT symmetry is precisely the class where the complex weight problem can be 
reduced to real weights, i.e., a sign problem. The spatial two-point functions of such models can 
exhibit three different behaviors: exponential decay, oscillatory decay, and periodic behavior. The 
latter two regions are associated with VT symmetry breaking, where a Hamiltonian or transfer 
matrix has complex conjugate pairs of eigenvalues. The transition to a spatially modulated phase 
is associated with VT symmetry breaking of the ground state, and is generically a first-order 
transition. In the region where VT symmetry is unbroken, the sign problem can always be solved in 
principle. Moreover, there are models with VT symmetry which can be simulated for all parameter 
values, including cases where VT symmetry is broken. 

PACS numbers: ll.30.Er, 11.10. Wx, 64.60.Bd, 64.60.De 
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I. THE COMPLEX WEIGHT PROBLEM 



The sign problem occurs in many branches of theoretical physics, including particle 
physics, many-body theory, statistical physics and condensed matter theory The prob- 
lem arises when the expected value of an observable is computed as a sum over non-positive 
weights, and perhaps should be called more generally the complex weight problem. For 
example, in lattice QCD, the introduction of a non-zero chemical potential makes the quark 
contribution to the functional determinant complex [Tj, [2] . Because the powerful simulation 
techniques of lattice gauge theory rely on a positive weight within the functional integral, 
progress in simulating QCD at finite density has been meager relative to what has been 
achieved in other aspects of QCD. Similar problems occur in other areas of physics, such as 
the Hubbard model away from half-filling [3] and systems with topological charges pQ. 

We will show that VT symmetry is a powerful tool for analyzing the sign problem. In par- 
ticular, "PT-symmetric systems are precisely the class of systems where the complex weight 
problem can be reduced to a sign problem, i.e., positive and negative weights. Furthermore, 
VT symmetry provides a natural classification scheme for the oscillatory behaviors observed 
in correlation functions in liquid and crystalline phases, giving us a simple, unified picture 
of behavior across many areas of physics. Finally, VT symmetry gives us important infor- 
mation about the solvability of the sign problem, and a framework for future progress in 
work on particular systems. 

The fundamental importance of VT symmetry was first pointed out by Bender and 
Boettcher in their seminal work on quantum- mechanical models [5]. Their work grew out of 
the observation that the Hamiltonian 

H=p 2 + igx 3 . (1) 

which arises naturally in the study of the Lee- Yang theory of phase transtions [6HTT]. has 
only real eigenvalues. Bender and Boettcher observed that the Hamiltonian H, while not 
Hermitian, is invariant under the simultaneous application of the symmetry operations parity 
V : x — > —x and time reversal T ■ % — > — i- This symmetry ensures that all eigenvalues 
of H are either real or part of a complex pair. The argument is simple: if H\ip) = E 
then HVT \ip) =VT H \ip) = VT E = E*VT Thus if E is an eigenvalue, E* is an 
eigenvalue as well. 
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II. QUANTUM MANY-BODY THEORY AND VT SYMMETRY 



We will now show that all quantum many-body problems involving a non-zero chemical 
potential may be described in terms of a non-Hermitian Hamiltonian with generalized VT 
symmetry. We start from a theory with a Hermitian Hamiltonian H and a conserved global 
quantum number N, obtained from a conserved current j^, that commutes with H. We 
assume that H is Hermitian and invariant under the combined action of time reversal T and 
a charge conjugation C that reverses the sign of j^. We take the number of spatial dimensions 
to be d — 1, and the spatial volume to be L d ~~ x . The grand canonical partition function at 
temperature T = and chemical potential // is given by Z = Tr[exp(—/3H + /3fj,N)]. 
If Z is written as a Euclidean path integral, the time component of the current j° will 
Wick rotate to ij d , while the chemical potential /i does not change. This leads directly to 
a non-positive weight in the path integral, and is the origin of the sign problem in finite 
density calculations: probabilistic methods do not work. The Euclidean space Lagrangian 
density may be written as C — i\xj d where C is the Euclidean Lagrangian for /x = 0; C — i\ij d 
is complex. The nature of the problem is changed by changing the direction of Euclidean 
time, so that we are now considering a problem at zero temperature with one compact spatial 
dimension of circumference (5. Upon returning to Minkowski space, j d does not rotate. We 
pick, say, the 1 direction to be the new time direction and take the new inverse temperature 
L to satisfy L ^> (3. When /i = 0, the original Hamiltonian is unchanged, but the partition 
function is now given by 



Tr 



(2) 



where 



Hp = H-ifi J d d - x x3 d . (3) 

The new Hamiltonian Hp is non-Hermitian, but possesses a generalized VT symmetry, where 
the role of V is played by the charge conjugation operator C that changes the sign of j° and 
N. Under the combined action of CT, j d — » —j d and i — > —i, leaving the Hamiltonian Hp? 
invariant. If we introduce the operator Hi = H — fiN, we have the relation 



Z = Tr 



Tr 



(4) 



under the space-time transformation that exchanges directions 1 and d. Note that Z is 
obtained from H L by a sum over all eigenstates, but is dominated by the ground state of 



Hp in the limit of large L. 

We have previously given an explicit example of the construction of Hp for the case of 
QCD with static quarks at finite density [12j. After the static quarks are integrated out of 
the functional integral, the effective action has the form 



Se ff = fd d x ^ 2 {F; U ) 2 - h F (e^Tr F (P) + e-^Tr F (P" 



(5) 



where P is the Polyakov loop P (x) = V exp i dtA^ (x, t) and hp is a known function of 
the temperature and the quark mass M, going to zero as M goes to infinity. In 1 + 1 di- 
mensions, the solution of this model reduces to solving a "PT-symmetric quantum mechanics 
problem on the gauge group. The Hamiltonian Hp, obtained from S e ff, is 

Hp = ^C 2 - h F [e^Tr F (P) + e"^Tr F (P + )] (6) 

and acts on the gauge-invariant physical states which are class functions of P obeying 
^ [P] = \& [gPg + \. The operator C 2 is the quadratic Casimir operator for the gauge group, 
the Laplace-Beltrami operator on the group manifold. The group characters form an or- 
thonormal basis on the physical Hilbert space: any gauge-invariant state \1/ can be expanded 
as \I/ [P] = Y^R a RTrR (P). The Hamiltonian H is not Hermitian when \x ^ 0, a direct man- 
ifestation of the sign problem, but it is VT symmetric when V is taken to be the charge 
conjugation operator C that changes to —A^. In figure [TJ we show the real part of 
the eigenvalues of Hp, measured in units where g 2 (3/2 is set to 1. The overall strength 
of the potential term is set by the dimensionless parameter 2h F /g 2 . In the upper graph, 
2h F /g 2 = 0.5, corresponding to antiperiodic boundary conditions for the heavy quarks. 
Note the formation of a complex conjugate pair of excited state eigenvalues as increases. 
It is also of interest to consider the case of periodic boundary conditions for the heavy 
quarks, corresponding to h F < [T3l - I15| . The lower graph shows the energy eigenvalues for 
2h F /g 2 = —0.5, where the ground state becomes part of a conjugate pair as /3/x increases. 
The coalescence of real energy eigenvalues into conjugate pairs with degenerate real parts is 
typical of PT-symmetric systems, and is usually referred to as broken VT symmetry. Note 
that this usage is somewhat different from what is meant by broken symmetry in Hermitian 
models, where only the symmetry of the ground state is considered. Note that for iV > 3, 
the heavy quark finite density problem of SU(N) gauge theory is in the universality class 
of the Lee- Yang problem for Z(N) spin systems [16]. 
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Figure 1: The real part of the SU(3) Hamiltonian Hp as a function of f3/i. The upper graph is 
for periodic boundary conditions for the heavy quarks, while the lower graph is for antiperiodic 
periodic boundary conditions. The energy has been scaled such that <? 2 /3/2 is set equal to 1. 

III. MODELS FROM STATISTICAL MECHANICS 



It is therefore not surprising that there is a class of PT-symmetric Z(N) spin models 
which are closely related. On each lattice site j there is a spin Wj, an element of the group 
Z(N) which may be parametrized as Wj = exp (27ririj/N) with rij G {0, 1, N — 1} denned 
modulo N so that and iV are identified. We take the operator V to be charge conjugation, 
acting as rij — > —rij, or equivalently Wj — > w*. The operator T is again complex conjugation. 
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Figure 2: Phase diagram for the d = 1 PT-symmetric Z(3) spin model in the Iir — hj plane at 
J = 0.2. The interpretation of regions la, lb, II and III are given in the text. 

Although V and T have the same effect on the w/s, one is a linear operator and the other 
antilinear. We will show below that V is implemented as a unitary matrix in the transfer 
matrix formalism. The Hamiltonian H is defined by 



w , 



(7) 



- pH = ^ - (wjiul + w*w k ) + ]T [h R (yij + wfj + hj ( 

(jk) j 

where /3 = 1/T, J, hi and hj are real and the sum over (jk) represents a sum over nearest- 
neighbor pairs. H is trivially PT-symmetric. This class of models has complex Boltzmann 
weights for N > 3 when hi ^ 0. 

We illustrate the rich behavior possible in these models using the case of a Z(3) model 
in d — 1. If hi = , then the transfer matrix T is Hermitian. When hi ^ 0, —/3H is no 
longer real and T is no longer Hermitian, but is VT symmetric. Figure [2] shows the phase 
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diagram in the Hr — hi plane for J = 0.2. There are four distinct regions. In region la, all 
three eigenvalues of the transfer matrix are real and positive. This region includes the line 
hi = 0, and has properties similar to those found in the Hermitian case. In region lb, all of 
the eigenvalues are real, but at least one of them is negative. In region II, the eigenvalue of T 
largest in magnitude is real, but the two other eigenvalues form a complex conjugate pair. In 
region III, the two eigenvalues largest in magnitude form a complex conjugate pair, and the 
third, smaller, eigenvalue is real. In both region II and region III, PTsymmetry is broken, 
but in different ways. Borrowing the terminology from PT-symmetric quantum mechanics, 
we will describe the behavior in region III as PT-symmetry breaking of the ground state, 
while region II is VT - symmetry breaking of an excited state. The behavior of the two-point 
function G (\j — k\) = (^w (j)w" f (k)^ differs substantially in the three regions. In region I, 
the two-point function falls off exponentially. We show typical behavior in region la in figure 
[3] for point A where (h^, hj) = (—0.45, 0.5). Similar behavior occurs in region lb, as shown 
in the figure for point B where (/ir, hj) = (—2.0, 1.5). Although the figure shows that the 
continuation of the two-point function away from integer values can be negative, note that 
the values at integer points are all non-negative. The two-point function at point C in region 
II where (h^, hi) = (0.25, 1.25) shows the damped oscillatory behavior associated with VT 
breaking in excited states. For the point D in region III, where (ha, hi) = (—0.5, 0.875), the 
VT breaking of the ground state leads to oscillatory behavior of the two-point function in 
the limit of large distance. Note that region III only occurs when ha is negative. For h^ < 
and hi = 0, the spin configurations with lowest energy have a two-fold degeneracy. With 
hi = 0, the ground state of the transfer matrix is unique. For the case hn < 0, hi = 0, and 
J large, the splitting of the two lowest eigenvalues of the transfer matrix in d = 1 is small. 
For sufficiently strong hi, the real parts of the two lowest eigenvalues of T merge, and VT 
symmetry breaking of the ground state occurs. 

VT symmetry is not always manifest in the Hamiltonian or transfer matrix. A simple 
criterion for a PT-symmetric Hamiltonian H (or transfer matrix T) has been given by Ben- 
der and Mannheim [T7]. If the characteristic polynomial det [H — XI] has real coefficients, 
then H has a generalized VT symmetry. In the case where H is real, interesting models 
arise when H is not symmetric. A notable example of this behavior is the ANNNI model 
|18| . The behavior seen above in the d = 1 Z(3) model is quite general. There are generi- 
cally three regions in "PT-symmetric models: region I, in which VT symmetry is unbroken; 
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Figure 3: The two-point function as a function of lattice spacing for the parameters corresponding 
to points A, B, C and D in figure [2] 

region II, in which VT symmetry is broken by one or more pairs of excited states becoming 
complex; and region III, in which VT symmetry is broken by the ground state becoming 
complex. In region II, thermodynamic properties are unaffected, but oscillatory behaviors 
appears in correlation functions. In region III, the system is in a spatially modulated phase. 
We wish to emphasize that the behavior of correlation functions seen in regions II and III 
cannot be obtained from classical spin models for which Hp is Hermitian: such behavior is 
incompatible with the spectral representation of the the correlation function for Hermitian 
theories. 



IV. ZEROS OF THE PARTITION FUNCTION IN REGION 3 

The change from region I to region II is generally not a phase transition in the thermo- 
dynamic limit. Within a transfer matrix framework, the largest eigenvalue is unique and 
real in both regions, so thermodynamic behavior is smooth. The change only appears in 
correlation functions. In condensed matter physics, the locus of points in parameter space 
where the change from region I to region II occurs has been known in statistical physics for 

8 



some time [19] and is often called the disorder line . The transition from region I or II to 
region III is generically a first-order phase transition. There is a general theory of parti- 
tion function zeros that can be applied to "PT-symmetric lattice models [20] • Under some 
technical conditions, the partition function in a periodic volume V = L d can be written as 

Z = £ e- pvf ™ + O (e- L / Lo e-P yf ) (8) 

m 

where / = min m Re [f m ] and Lq is of the order of the largest correlation length of the system. 
The / m 's have the interpretation of complex free energy densities, and are independent of L. 
These phases are stable if Re (f m ) — f or metastable otherwise. The zeros of the partition 
function are within O (e~ L ^ L °) of the solutions of the equations. 



Re(f m ) = ite (/„) = / 

Im(f m ) = j m (/ n ) + (2p + l)~ 

for some m ^ n and p G Z . We can apply this directly to region III, using the representation 

Z = J2 e- LEr + £ (e- LE * + e~ LE i) (9) 

r p 

of the partition function. We identify LE and LEq as f3L d f and f3L d fQ, so that the 
partition function has a zero for values of the parameters such that 

flMAH^^ do) 

This tells us that the zeros of the partition function lie on the boundary of region III, defined 
by Im [/o] = 0, in the limit V — > oo. As the volume of the system is taken to infinity, the 
zeros of the partition function lie asymptotically on the boundary between phases. Note 
that this analysis depends on L remaining finite. At a 2nd-order transition, L goes to 
infinity and the approximation is invalid. 



V. TOWARDS A SOLUTION OF THE SIGN PROBLEM 

The difficulty presented by the sign problem depends directly on VT symmetry breaking 
or its absence. Mostafazadeh has proven that when VT symmetry is unbroken (region I), 
there is a similarity transformation S that transforms a PT-symmetric Hamiltonian H into 
an isospectral Hermitian Hamiltonian via = SHS^ 1 [21J. This eliminates the sign 



9 



problem for "PT-symmetric quantum Hamiltonians throughout region I, if S can be found. 
This theorem also applies to "PT-symmetric transfer matrices, but a further restriction to 
positive eigenvalues is necessary for the elimination of the sign problem. 

In regions II and III, the sign problem has an underlying physical basis. The negative 
weight contributions to the partition function Z arise from the contributions of complex 
conjugate eigenvalue pairs associated with VT symmetry breaking. It is that breaking that 
in turn gives rise to the oscillatory and damped oscillatory behavior of two-point functions 
characteristic of many physical systems. However, there are VT symmetric models like the 
ANNNI model [18] where the classical Hamiltonian, corresponding to the action in the path 
integral formalism, is real. Such a model can be simulated with no difficulties of principle 
throughout its parameter space. The key point is that the the transfer matrix is real but 
not symmetric. The existence of an antiunitary involution commuting with the Hamiltonian 
implies that there is a basis in which H is real; see, e.g., j22]. This theorem easily extends 
to the case of those PT-symmetric systems for which {VT) 2 = 1, and can be applied to 
transfer matrices as well as Hamiltonians. This suggests the existence of a broad class of 
PT-symmetric models which can be simulated in all three regions, but the extent of this 
class is as yet unknown. 
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